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Abstract 

The (2+l)-dimensional Thirring model is studied by using the Gaussian ap- 
proximation method in the functional Schrodinger picture. Although the dy- 
namical symmetry breaking does not occur in the large N limit, it does occur 
in the Gaussian approximation which includes the higher order contributions 
in 1/N. 
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I. INTRODUCTION 



The Gaussian variational method in the functional Schrodinger-picture has been shown 
to be useful for the study of quantum structures of field theories . This method has been 
shown to be well-suited not only for the study of bosonic field theories, but also for the study 
of fermionic field theories |@] . The Gaussian method is non-perturbative from the viewpoint 
of both the ordinary weak coupling expansion and also the 1/N expansion, and hence it can 
give the better informations than the large-N expansion in principle, although it has not 
been realized yet. It is the purpose of this paper to give an example where the Gaussian 
approximation method provides better information than the large N approximation. 

In the Floreanini- Jackiw representation the fermion field operators , which satisfy the 
anticommutation relations 



are realized as 



{ipi (x, t) , V^] iy,t)} = iSijS{x, y), (1) 



= —=[u{x) + 



V2' ' ' 5nt(x)^' 

with anticommuting Grassmann variables u and u"''. For the variational approximation, we 
take the trial wavefunctional in the Gaussian form 

I ^ >^\G>=j^-^^^exp[J^^u^x) G{x,y) u{y)], (3) 

< ^ |^< G 1= ^^^1 explj^ u\x) G{x, y) u{y)l (4) 



where 



G = {G^y\ (5) 

This prescription has been used in to study the Gross-Neveu model. In the present 
paper, we apply this Gaussian variational method to study some aspects of the (2 + 1)- 
dimensional Thirring model. In the models considered so far, the Gaussian variational 



method gives qualitatively the same results as the large approximation. We will show 
that the Gaussian approximation method applied to the (2+l)-dimensional Thirring model 
provides new informations that can not be obtained in the large limit. 

The massless Thirring model has, at the classical level, parity and chiral symmetry in 
two and four component representations of fermions, respectively. Recently it has been 
shown [^,^ that in this model, dynamical breaking of parity (chiral) symmetry occurs as 
a cooperative effect of different orders in expansion. In this paper we consider this 
problem in the context of Gaussian approximation scheme. It is shown that the symmetry 
breaking occurs in the Gaussian approximation although it does not in the large N limit. 
This agrees quahtatively with the results of @,|^. 

II. DYNAMICAL SYMMETRY BREAKING IN TWO-COMPONENT 

REPRESENTATION 

A. Schrodinger picture Gaussian approximation 

The (2+l)-dimensional massless Thirring model is described by the Lagrangian density 

where a = 1, 2, A^. In (2+1) dimensions the 7 matrices for two-component representations 
of fermions can be represented by the 2x2 Pauli matrices, 

/ = a3, ^'=za\ ^^ = ia\ (7) 

As is well-known, the model has parity invariance at the classical level, under which the 
fermions transform as 

?/'(x°, x^,x^) 7V(a;°, -x\x^). (8) 
The mass term ipip changes sign under the transformation and thus breaks the symmetry. 



3 



In order to compute the effective potential and see if dynamical symmetry breaking 
occurs, we introduce the vacuum condensate of fermion bilinear, 

(^a = -^< 'ipaA >, (9) 

where a is a color index, and compute the effective potential as a function of aa- One can 
achieve this by writing the Lagrangian as 

C = Co + aa{cTa + ^Ipa'^a)- (10) 

By using the Gaussian trial wave functional @ and (1), we calculate the vacuum expec- 
tation value of the Hamiltonian of the model, 

<G\H\G> = ^ Jd\d'yTr[i-t^'Yd.ix,y) + |/)7 V + |«a(x)5(x, y)7°)^^(l/, x)] 

- ^Tr[^'Yn{x, yh'iMy, x)6{x, y)] 

+ J (i^{^Tr[7V^(^,a;)]rr[7Vfi(x,x)] (11) 



where 



n{x,y) =2 <G\^^^G>=<x\{l + G)S-\l + G)\y>, (12) 
S=G + G, (13) 
d,{x,y)^^{6{x^-y')), (14) 

and Tr denotes the trace over Dirac spinor and color indices. In order to calculate the 
Gaussian approximation, it is convenient to introduce the current = < ip'^f^ip >, 
which can be realized by introducing the Lagrange's multiplier into Hamiltonian, 



Heff = \j d'^xd^yTr[h^{x,y)Q{y,x)] 



+ I d'x{^^A^{x)A,{x) + (3,{A^{x) - -^Tr[^'YQ{x, x)]) + a.aj (15) 



where 



hn{x,y) = -i-f^-f'di{x,y) 



+ ^7°7'^(2/(x, y) - n{x, y))5ix, yh\ + |«a(x)<5(x, yh'. (16) 
The effective potential is obtained by minimizing H^ff with respect to the kernel G{x,y). 



B. Effective Potential in the Large N Limit 

Since the second term in the hn of Eq(|16|) is of order, the effective potential in the 
large N limit is obtained by considering 

H,ff = ^Trh^n + Yg^^^i^ + + «acra, (17) 

where 

^ -ii'i'di{x, y) + y)7° - -^M^)i'r5{x, y). (18) 



In order to obtain the effective potential of the system, we take variations on < G \H\ G > 
with respect to G and G . The invariance oi < G \H\ G > under these variations yields the 
equations 

[I - G)hN{I + G) = Q (19) 
(/ + G)/i^(/-G) =0. (20) 

These two conditions are shown to be equivalent under the condition of (^, and can be 
solved by the well-known method 0. The Eq.(|l^) can be rewritten in the form 

hl-Kl+[h,KN] = Q (21) 

where and Kn = IinGn are 2x2 matrices in the Dirac spinor space. Any 2x2 matrix 
can be expressed linear combination of r*^: 

ro = j, T' = -t^\ r' = -t^', r=^ = 7°. (22) 

K]\f and is then decomposed as 

Kn{x, 2/) = E / 7|Tie-^^-(^-^)ir;va(p) (23) 



a=0 

72, 



- {P2 + gP2)r' + {vi + ^7/5i)r' + ga^T^ (24) 



The Eqs.(||),(||) and (|2T]) yield 



KnqKni = gPoiP2 + 9^2)+ i{Pi + gPi)KN3 - igaaKN2, 

KnqKn2 = -gPoiPi + 9Pi) + i{P2 + gp2)Km + igaaKm, (25) 



KmKm = -g'^Poota - i{p2 + gP2)KN2 - i{pi + gPi)K^ 



Nil 



where the summation convention for the index i(=l,2,3) is imphed. The Eq. (^5D has the 
non-trivial solutions 



Kno = ±^/hihi, 

KNi = hohi/KNo (26) 
as well as the trivial solution Km = ±gPo, and this gives the solution for the kernel 

GN{x,y) = {h]^^KN){x,y) 

Following the similar procedure as in the case of Gross-Neveu model 0, the effective 
potential K//, which is defined by 

iJe// = / dWejj, (28) 

can be determined as 



72 

K// = N jd-'xi-j -7^,{gPo + ^-{pi + gPi?-{P2 + gP2? + g^a^} 



+^d^A^ + aa] (29) 



2Ng N' 



2tt 2tt Arc Qn 

+^A''A^ + /3>'A^ + Naa, (30) 
2g 



where 
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Ji = J pdp 

h = Jp^ dp. (31) 

Here we have taken the solution with minus sign in Eq.(^) and (pTf), which corresponds 
to lowest energy. Since we want to see whether dynamical symmetry breaking occurs, the 
auxiliary field variables in Eq.(^) should be eliminated via variational method. If we 
neglect the irrelevant infinity, the effective potential becomes 

Veff = -^gWlo{M) + ^g'\a f +Naa. (32) 

47r DTT 

By eliminating a, the effective potential can be written as 

N 4 
V,ff = --9'a\lo{M)--g\a\). (33) 

This effective potential can easily be renormalized in the form, 

Veff = -^9>^ (34) 

where = g'^Io{M) is the renormalized coupling constant. This implies that the symmetry 
breaking does not occur in the large N limit. 

C. Beyond large N limit 

To compute the full Gaussian effective potential, it is convenient to rewrite the original 
effective hamiltonian density (|TBp in the form 

H^ff = ^Tr[hn] + -|-Tr[7^''^]7^^^]] + ^^A'^A, + /3,A^ + aaa^ (35) 



where 



h^hr, + ^7 V(/ - ^{x, x)h\. (36) 



Taking variations on H^ff with respect to G and G yield 



(J - G)h{I + G) = 0. (37) 

This equation has the same form as Eq.(^) with replaced by h. Therefore we can solve 

h^-K^ + [h, K] = (38) 

where K = hG. We can follow the same procedure as in the case of large N-limit and obtain 
the nontrivial solution for K 

K{p) = -^^{h.h^T^ + hoh.P), (39) 



and 



1 

Vhihi 



G{p) = -^==h,P, (40) 



which corresponds to the minimum energy solution. We also obtain the solution for Q, 
Q = I + G, from the relations among G, G\ K, K\ h, and /i^. Then using the Eq. (|36|) one 
finds 



where 



G^iO) - / ^Mp)- (42) 



From Eq.(iy) and (^T]), we obtain the consistency condition 



and therefore, 



Gi(0)=G2(0) = 0, 

n fn\ f 3 Pi ^2 , / 9 02 ^2 , 21-i 
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-TT-ih-m), (44) 



where 

m^^-a-^GM- (45) 
Extremizing the effective potential with respect to 6*3(0), we obtain 

Veff = -— /o M m2 + — m 3 +J^iVV + — 
47r 67r 12 S^f 

—am + A^acr (46) 

Eg. (^51) can then be written as 

2777, 

a - — + 3m27rA^(/o - m) = 0, (47) 



which can be understood as ^K// = 0. 



By ehminating a by ^K// = and a by ^K// = 0, the effective potential becomes 

=^{N~ ^g)m' + —\m\\ (48) 

which needs to be renormalized. Vgff can be made finite by defining the renormalized 
coupling constant Qr by 

-l^V-^,). (49) 
9r 9 47r 

Then, the renormalized effective potential cab be written as 

V,ff = -^m' + f\m\'. (50) 

6gr OTT 

This implies that the symmetry breaking occurs ii Qr > 0- Since a = Nm/3 in the extrimum, 
the second term of Eg . (|50|) is smaller than the first term by Eq.(^) shows, in a 

straightforward manner, why symmetry breaking phenomenon cannot be seen in the large 
N limit. 
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III. DYNAMICAL SYMMETRY BREAKING IN FOUR-COMPONENT 

REPRESENTATION 



To restore parity symmetry even with the exphcit mass term, we consider four-component 
representation of Dirac spinor. This can be achieved by considering a doublet of two com- 
ponent spinors, 



(51) 



for which we define the 4x4 7-matrices 



=iT 



13 = « 



-(Ti 



2 -n -n 5 -p 

7 = «r23, 7 = «r33, 7 75 = ro2. 



(52) 



where Fjj = Oi® aj. In this representation, mass term does not spoil parity symmetry, but 
break the continuous U{1) x f/(l) chiral symmetry 0, 



(53) 



down to the f/(l) subgroup. 

Thus the existence of mass term signals the breaking of such chiral symmetry. We 
adopt the same procedure, as in two-component case, to obtain the effective potential, and 
introduce two vacuum condensates of fermion bilinear 



g — ig — 

a = -- Ki/jijj >, TT = — < V75V' > • 

Then the effective Hamiltonian can be written as 



(54) 



(55) 



where 



h 



-i-i^l'di{x, y) + ^a{x)5{x, t/)7° • 

/v iV 
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+ ^7V(1 - ^^(x,x))7%, - |r7%5. 



(56) 
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By following the same procedure as in the two component case, we obtain the effective 
potential 

where 0^ = cr^ + vr^. The model exhibits the dynamical breaking of the chiral symmetry in 
the exactly same manner as in the two-component case. 



IV. CONCLUSION 

The Gaussian approximation method has been shown to be useful in studying the non- 
perturbative aspects of quantum field theories. In the models considered so far p|, the 
Gaussian approximation method gives qualitatively the same results as those of large N 
approximation. This is partly due to the fact that the main non-perturbative feature of 
those models appears at the leading order in expansion. On the other hand (2+1)- 
dimensional Thirring model has the dynamical symmetry breaking not in the leading order 
of N,but in the combination of different orders in 1/A^ 0)0- In two component formalism, 
this non-perturbative phenomena, from the point of view of large N expansion, occur in 
the range of qr, < Qr < gc = TQ^^Pi~^l^) HI- this paper we use the Gaussian 
approximation scheme to reveal this. We have shown that in this scheme the dynamical 
symmetry breaking does occur as a cooperative effect of the leading and next-to-leading 
orders in 1/A^. Since the Gaussian approximation gives better results than the large N 
approximation, we believe that the Gaussian approximation is well-suited for the difficult 
non-perturbative problems such as the existence of bound state. 
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